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Abstract 



By virtue of advances in next generation sequencing technologies, we have access to new genome sequences almost daily. The tempo of these 
advances is accelerating, promising greater depth and breadth. In light of these extraordinary advances, the need for fast, parallel methods to 
define gene function becomes ever more important. Collections of genome-wide deletion mutants in yeasts and E. coli have served as 
workhorses for functional characterization of gene function, but this approach is not scalable, current gene-deletion approaches require each of 
the thousands of genes that comprise a genome to be deleted and verified. Only after this work is complete can we pursue high-throughput 
phenotyping. Over the past decade, our laboratory has refined a portfolio of competitive, miniaturized, high-throughput genome-wide assays that 
can be performed in parallel. This parallelization is possible because of the inclusion of DNA 'tags', or 'barcodes,' into each mutant, with the 
barcode serving as a proxy for the mutation and one can measure the barcode abundance to assess mutant fitness. In this study, we seek to fill 
the gap between DNA sequence and barcoded mutant collections. To accomplish this we introduce a combined transposon disruption-barcoding 
approach that opens up parallel barcode assays to newly sequenced, but poorly characterized microbes. To illustrate this approach we present a 
new Candida albicans barcoded disruption collection and describe how both microarray-based and next generation sequencing-based platforms 
can be used to collect 10,000 - 1,000,000 gene-gene and drug-gene interactions in a single experiment. 



Video Link 



The video component of this article can be found at http://www.jove. com/details. php?id=2864 



Protocol 



1. Background information 

There are several ways to generate mutants that carry barcode tags. The current gold standard is the Yeast KnockOut (YKO) collection created 
by a consortium of labs and completed in 2002 1 . Since the original YKO was introduced, other yeast collections have been generated; in different 
strain backgrounds, using over-expression constructs, and in other microbes such as E. coli 2 . In parallel, the effort to create barcoded shRNA 
libraries is proceeding rapidly, and in fact, many of the design principles for these mammalian collections have been adopted from yeast. To 
demonstrate how barcoded transposons can be a rapid, widely applicable strategy for creating systematic mutant collections, we focus on one 
collection we recently created in the human fungal pathogen, Candida albicans. Our work on Candida was based on the success of barcode 
screens in S. cerevisiae, and is used here as an example organism. The sample protocol, can with minor modifications be used to screen any 
organism that can be grown in suspension culture. Because few organisms have the requisite high rates of transformation and efficient mitotic 
recombination needed to create perfect deletion mutants, accordingly we developed a protocol that uses transposon mutagenesis in vitro to 
mutagenize a genomic DNA library, and then transformed these barcoded genomic fragments into Candida albicans 3 - 4 . Inspired by the success 
of the original YKO collection and its role in fundamental discoveries on the nature of gene networks 5 ~ 8 , genome-wide haploinsufficiency 9 , drug 
target and mechanism of action 10 ' 11 , and the essentiality of all genes in the genome 12 we anticipate broadening this approach to other microbes 
will be extremely fruitful. 

The protocol below assumes that the desired mutant collection has been created (e.g. YKO or Candida albicans disruption collection) and is 
available as individually archived strains. For a detailed description of strain construction see 1 ' 13 U . 

2. Combine individual mutants into a single pool 
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1 . Allow one week to generate pooled cell aliquots (can be stored indefinitely at -80°C). 

2. Thaw the frozen glycerol stocks for the strains of interest completely but do not let cells remain thawed for >2hrs. 

3. Sterilize a 96-well pin tool, dip the pin tool in water to remove any remaining cells, followed by 2 dips in 70% ethanol baths (e.g. pipette tip 
box lids), flame the pin tool and cool for 1 minute. Take care to flame the pin tool away from the ethanol baths. The level of the ethanol baths 
should exceed the level in the water bath to ensure all carry-over cells are flamed and removed. Replace water every 4-6 pinnings. 

4. Insert the sterile 96-well pin tool into a thawed 96-well plate, swirl gently and transfer cells to a Nunc Omni Tray containing YPD-agar 
including the appropriate antibiotic. Grow colonies until they reach maximal size at 30°C (2-3d). To conserve plates, we find it most useful to 
consolidate four 96 well plates onto a single Omni-tray with "384 strains. 

5. After colonies have grown, note any missing or slow growing strains and repin these at "2X the cell mass as the rest of the strains. 

6. Working in a microbiology environment (with flame and sterile labware, flood plate with 5-10 ml media, soak for 5 min and resuspend colonies 
with a cell spreader. Pour the liquid plus cells into a 50 ml conical centrifuge tube, and add glycerol to 15% or DMSO to 7% (vol/vol). 

7. Measure the OD600 of the pool and adjust (by dilution or centrifugation) to a final concentration of 50 OD6oo/ml with media containing 1 5% 
glycerol or 7% DMSO. 

8. Aliquot in 40 pi volumes in PCR strip tubes and freeze at -80°C. 

3. Experimental pool growth 

This procedure is outlined in Figure 1. 

1 . Thaw pool aliquots (in PCR tubes) on ice. If not using robotics, skip to step 5. 

2. Immediately dilute (gently!) the pool into media with drug or condition of choice, inoculating at an OD600 of 0.0625 in a total volume of 700 pi 
in a 48-well plate. Include at least one appropriate solvent control on the plate. For experiments that extend beyond 5 generations of growth 
(i.e. more than 1 well), fill adjacent wells with media or condition of choice, but NO cells. 

3. Seal with a plastic plate seal; if the condition requires aerobic growth (e.g., non-fermentable carbon sources), use a 21 gauge needle to 
pierce the holes in the membrane seal towards the side of each well. 

4. Grow in a spectrophotometer, shaking at 30°C with an experimentally determined shaking regimen (e.g. shake 14 min at highest setting (or 
temperature-controlled shaker), read wells, resume shaking). Part of the cell suspension can be harvested by the robot and saved on a cold 
plate on the robotic deck at user-defined generation time points, typically 5, 10, 15 and 20 generations of growth. 
(http://med.stanford.edu/sgtc/technology/access.html, for details contact C. Nislow or G. Giaever). 

5. For manual cell growth, inoculate a 50 ml culture at a starting OD600 of 0.002 in a 250 ml culture flask. Shake at 30°C at 250 rpm until cells 
reach a final OD600 of 2.0 (for Saccharomyces or Candida) for "10 generations of growth. Additional generations of growth can be obtained by 
diluting cells at an OD600 of 2.0 back to 0.02 in a fresh flask. 

6. Harvest at least "2 OD600 units of cells for each sample/time-point into Safe-Lock microcentrifuge tubes. 

Note: Always collect a starting cell sample (i.e. a "TO time point") to assess initial strain representation in any newly created pool by adding 1-2 
OD600 of pool directly from the freezer aliquots to a 1.5 ml tube and processing as described below. 

4. Genomic DNA extraction, PCR and microarray hybridization or sequencing 

1 . Purify genomic DNA from "2 ODeoo of cells with the Zymo Research YeaStar kit according to manufacturer's instructions (Protocol I if using 
yeast DNA), or another suitable method specific to the organism of interest (standard phenol/chloroform extraction followed by alcohol 
precipitation works well for diverse microbes). If using the YeaStar kit, elute the DNA with 300 pL of 0.1X TE instead of the 60 pL of 1X TE 
specified in the protocol. Genomic DNA can be stored indefinitely at -80°C. 

2. Set up two PCR reactions for each sample, one for the uptags and one for the downtags, with the reaction conditions as follows: 33 pi ddH20, 
6 pi 1 0X PCR buffer without MgCI 2 , 3 pi 50 mM MgCI 2 , 1.2 pi 10 mM dNTPs, 1.2 pi 50 pM Up or Down primer mix, 0.6 pi 5 U/pl Taq 
polymerase, "0.1 pg genomic DNA in 15 pi. Total volume is 60 pi. Thermocycle under the following conditions: 94°C 3 min, 30 cycles of 94°C 
30s, 55°C 30s, 72°C 30s; then 72°C 3min, and hold at 4°C. Check the resulting PCR products on a gel; a 60 bp product for both PCRs is 
expected for amplicons used for hybridization and 130bp for amplicons for barcode sequencing). The PCR products can then be stored at 
-80°C indefinitely. 

3. Prewarm the hybridization oven temperature to 42°C and set up a boiling water bath and ice bucket containing an ice-water slurry. 

4. Pre-wet the arrays by slowly filling with 120 pi 1X hybridization buffer. 

5. Incubate in the hybridization buffer at 42°C and 20 rpm for 10 minutes. 

6. Prepare 90 pi of hybridization mix per sample, plus one extra as a buffer, as follows: 75 pi 2X hybridization buffer, 0.5 pi B213 control 
oligonucleotide (0.2 fm/pl), 12 pi mixed oligonucleotides (12.5 pm/pl), 3 pi 50X Denhardt's solution) in lock-top 0.5 ml tubes. 

7. Add 30 pi uptag PCR and 30 pi downtag PCR to 120 pi hybridization mix for a total volume of 150 pL. Boil for 2 minutes and set in ice-water 
for at least 2 minutes. Briefly spin the tubes prior to use. 

8. Remove the pre-hybridization buffer from the arrays and add 90 pi hybridization/PCR mix. To prevent evaporation, cover the array gaskets 
with a Tough-Spot. Hybridize for 16 hours at 42°C, 20 rpm. 

9. Freshly prepare 600 pi biotin labeling mix per sample plus one extra, as follows: 180 pi 20X SSPE, 12 pi 50X Denhardt's, 6 pi 1% Tween 20 
(vol/vol), 1 pi 1 mg/ml streptavidin-phycoerythrin, 401 pi ddH20. Store all streptavidin-PE samples in the dark. Aliquot 600 pi into 2 ml tubes. 
Remove Tough-Spots from chips. 

10. Slowly remove hybridization mix from the arrays by pipette and fill microarrays with 120 pi Wash A. Prime the Affymetrix fluidics station. 

11 . Wash the arrays using an Affymetrix fluidics station according to the manufacturer's instructions, using "Gene-Flex_Sv3_450" protocol with 
the following modifications: 1 extra step with Wash A (1 cycle, 2 mixes) before staining, Wash B temperature 42°C instead of 40°C, stain at 
42°C instead of 25°C. It is also possible to perform the post hybridization wash, the biotin staining, and the post staining wash manually (see 
p. 396 in reference 15). Following fluidics operations, run the fluidics station "SHUTDOWN_450" protocol. 

12. Following the wash, ensure no air bubbles are present. If necessary add 90 pi Wash A and pipette slowly until the bubbles disappear. If there 
are any marks or smudges on the array surface, clean the glass window with isopropanol and a lint-free tissue. Apply fresh Tough-Spots over 
the gaskets/septa and place the arrays in the scanner. 

13. Scan in an Affymetrix GeneArray scanner at an emission wavelength of 560 nm. 

5. Array analysis (see Figure 2 for an example obtained using the Candida albicans 
disruption collection) 
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1 . Outlier masking: Because the Affymetrix TAG4 array contains 5 replicates of each barcode complement dispersed randomly any barcode 
probe that deviates from the replicates can be masked and discarded. To do this, for each array feature that appears to be an outlier based 
on its signal compared to the 4 other replicates of this feature, our software first examines the 5 features that surround the suspect outlier, 
producing a matrix of 25 features with the suspect feature at the center. If >1 3/25 probes in this region differ from each of their individual their 
trimmed replicate mean (the mean of the middle three replicates, excluding the highest and lowest replicates) by more than 10%, this probe is 
then discarded from further analysis. Because such outliers are most often the result of post-hybridization washing inconsistencies, we 
expand or "pad" the region containing suspect probes. Pad such probes by including all probes within a 5-probe radius, as defined by ((X1-X2) 2 
+ (yi-y2) 2 ) ,/2 < 6 where xi, X2, yi, and 72 are the x and y coordinates for the two features. Finally, discard features for which standard deviation 
(included in the .eel file for Affymetrix arrays) of feature pixels/mean feature pixels. After removing outliers, average the intensity values for all 
remaining replicates. 

2. Removing unusable tags: Tags with low-intensity values will give poor-quality results and must be removed. An exclusion cutoff for these low 
intensity probes can be calculated as follows: 

a. For any treatment-control pair of arrays, calculate \og2((h-b g )l(irb g )) for each tag, where i c is the control intensity, it is the treatment 
intensity, and b g is the mean intensity of the unassigned tag probes. 

b. Pair the uptag and downtag ratios by strain and for each tag pair, take the minimum intensity for the two tags in the two samples. Sort 
the ratio pairs by this minimum intensity. 

c. Use a sliding window size of 50 on the ranked ratio pairs, calculate the correlation of uptag and downtag ratio pairs within the window. 
Also calculate the average of the minimum intensities calculated in the previous step. 

d. Slide the window by 25 pairs, and repeat the previous step until all pairs have been crossed. 

e. Plot the average minimum intensity versus the uptag-downtag correlation for all windows. 

f. Finally, choose an intensity threshold; generally we use the intensity value where the correlation first reaches 80% of its maximum 
level. Flag and remove from further analysis any tags below this cutoff. 

3. Saturation correction: Because each feature on the barcode microarray can become saturated, the signal on the TAG4 array is not linearly 
related to tag concentration. To correct for this saturation follow the protocol described in reference 16 . 

4. Array normalization: For each array, normalize the uptags and downtags separately. To quantile normalize, rank the values obtained from 
each array for uptags and downtags in order of increasing intensity. To mean normalize for each set of uptags and downtags, divide by the 
mean. The mean normalization of all arrays is a second step following transforming the raw data to the mean of each array. 

5. Calculating sensitivity scores for control-treatment comparisons: To use log2 ratios as a metric of sensitivity: For each strain, calculate log2((/J c 
-bg)/(fj t -bg)), where u c is the mean intensity for the control samples, u t is the mean intensity for the treatment samples, and b g is the mean 
intensity of the unassigned probes. Strains with a positive log2 ratio are sensitive to the treatment, and those that are resistant have negative 
log2 ratios. 

6. Assessing fitness of barcoded yeast strains by sequencing 

1 . Isolate DNA from the deletion pools as described for microarrays. 

2. Amplify each 20-mer uptag barcode with composite primers comprised of the sequences of the common barcode primers and the sequences 
required for hybridization to the lllumina flowcell (See Table of lllumina primers and Figure 3 for a diagram of the amplicon). These primers 
can be used desalted without additional purification. PCR is performed in 100uL, using Invitrogen Platinum PCR Supermix (Cat. No. 
11306-016) with the following conditions: 95°C/3 min; 25 cycles of 94°C/30 sec, 55°C/30 sec, 68°C/30 sec; followed by 68°C/10 min. 

3. Purify the PCR product (~130bp) with Qiagen MinElute 96 UF PCR Purification Kit (Cat. No. 28051). 

4. Following PCR purification, quantify DNA with the Invitrogen Quant-iT dsDNA BR Assay Kit (Cat No. Q32853). Do not rely on 260/280 
readings! 

5. Normalize DNA concentration to 10ug/ml and pool equal volumes of normalized DNA. 

6. Separate pooled DNA on a 12% polyacrylamide TBE gel for 3-4 hours depending on voltage used. 

7. Stain gels with ethidium bromide (Sybr Green should work as well) for 30 minutes. 

8. Locate the band of interest on a longwave UV lightbox (wearing appropriate face protection), cut it out and extracted the DNA using the crush 
and soak method 17 followed by ethanol precipitation. 

9. Confirm that the appropriate size DNA (130bp) has been isolated and that primers have been removed using the Agilent Bioanalyzer High 
Sensitivity DNA kit (Cat No. 5067-4626). 

10. Sample sequencing: 

a. lllumina GAIIx platform: 

1. Generate clusters on a Single-Read flowcell using the cBOT and Single-Read Cluster Generation Kit (Cat No. GD-300-1001). 
For Read 1, UP and DOWN-tag modified sequencing primers are pooled at a 100uM stock concentration and added to a 
strip-tube (0.6uL of each 100uM sequencing primer in 120 pi HT1). Recipe SR_Amp_Block_StripTubeHyb_v7.0 is used to 
generate R1 clusters. 

2. Sequencing on the Genome Analyzer llx. Following 18 sequencing cycles, the Paired-end module is used to strip the 
synthesized first strand and rehybridize the flowcell, using the lllumina R1 (below). Clusters are regenerated and sequenced for 5 
cycles to capture the index tag. 

11. The Index tag sequence is used to bin sequences into experimental bins. 

12. Within each experimental bin, the yeast barcode sequences are tallied to give a total number of counts for each barcode. 

13. Counts are quantile normalized so each experiment has the same count distribution. By analogy with barcode microarray fitness experiments, 
fitness defect ratios for each strain are calculated and expressed as the log2 ratio (control/treatment). Positive fitness defect scores signify a 
decrease in strain abundance during drug treatment and suggest that the wild-type version of the gene deleted in that strain is required for 
resistance to that drug or inhibitor. 

Note: Considering Bar-seq as an alternative to array hybridization. With the costs of high-throughput sequencing decreasing, using 
high-throughput sequencing as a readout of tag abundance is becoming feasible and in many cases, is more cost effective 18 . In this way, 
amplified PCR product is measured directly as 'counts' rather than as signal intensity as hybridized to an array. This eliminates false negatives 
and positives that arise from tag cross-contamination, saturation, or issues arising from very high or very low signal intensities. Furthermore, 
multiple experiments can be combined prior to sequencing by the addition of a 4-8 base DNA index 19 . Because the yeast barcodes are 20 bp, 
one single, 2-step read of 26-28 bases captures both multiplex index and unique barcode , allowing for extreme 1 00+ multiplexing. At the time of 
this writing, Bar-seq offers a cost advantage over Bar-code microarrays, and furthermore, Bar-seq is inherently flexible such that as the number of 
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reads/run increases, the level of multiplexing can increase to further decrease costs. Several "mid-capacity" sequencers from all of the major 
platform manufacturers will further democratize Bar-seq, with sequencing likely to become the readout of choice. 

This protocol has also been validated on the lllumina HiSeq2000. 

An excellent demonstration of the use of Bar-seq to address a fundamental biological question in Saccharomyces cerevisiae growth control is 
presented in a recent study by Gresham et al. 20 who outline several important experimental design and interpretation guidelines. 



7. Validation of pooled screening data 



The results from any functional genomics screen should be verified using the individual strains in isolated culture. Because each experiment will 
differ in terms of the number of sensitive strains, selecting the number of candidate strains to confirm is somewhat arbitrary. As a guide, ranking 
the most sensitive strains by the log2 ratio or z-score and testing the top 25-50% of the candidates (which typically translates to 2-3 standard 
deviations from the mean of all strains in the pool) is a good balance between costs and benefits. Individual confirmations can be performed in 
any flask but we perform these tests for 5 generations of growth in 96 well plates using a starting inoculum of 0.06 OD600 in 1 00 pi of media in a 
shaking spectrophotometer, taking measurements every 15 minutes (See Figure 4). 



8. Representative Results 



Once a genome-wide screen is complete, and the arrays have been normalized and the behavior of each strain compared to a control treatment 
(e.g. by comparing microarray intensities or sequencing counts/strain) the data are most easily manipulated in an excel file with genes ranked by 
the log2 ratios of control/experiment. In this manner, the greater the negative log2 ratio, the more sensitive that particular strain is to the test 
condition. These excel files can be plotted in a variety of graphing software packages. We find it simplest to plot the log2 ratios on the Y axis and 
the gene or ORF names on the X axis. In the example shown in Figure 2a, such a plot of clotrimazole treatment (a known antifungal agent) is 
shown. All strain that are significantly sensitive to treatment with a log2 ratio of 2 are highlighted in red, and we would typically verify many of such 
strains in individual growth assays of each mutant in the presence of the same concentration of drug. In this example, 4 strains are highlighted, 
NCP1 , ERG2 and 2 independent alleles of ERG11, the known protein target of clotrimazole. Each of these 4 genes is directly involved in 
ergosterol biosynthesis, the yeast equivalent of cholesterol. For example, NCP1 encodes a NADP-cytochrome P450 reductase that is involved in 
ergosterol biosynthesis and which is associated and coordinately regulated with Erg11 . This example highlights the fact that the known drug 
target (Erg11 ) is identified in this unbiased screen, as well as several other key components of the target pathway. Finally, several of the 
highlighted genes in red represent genes that may be involved in ergosterol biosynthesis or in distinct biological processes. As mentioned above, 
each strain detected as sensitive in a pooled screen should be verified as to its sensitivity in individual growth assay. In the example shown in 
Figure 2b, four strains are confirmed to be sensitive to clotrimazole based on their decreased growth relative to the wild-type parent strain, 
BWP17. These individual growth curves highlight an important feature of such pooled gene-drug screens; that is the absolute rank of a particular 
strain does not necessarily reflect its exact level of sensitivity. Furthermore, Figure 2b also shows the value of having multiple alleles for each 
gene, in this case, the two erg11 disruption mutants have slightly differing sensitivities. Correlating the nature of these disruptions with the degree 
of sensitivity can provide additional insight into the drugs mechanism of action. 
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using an aliquot of 
pooled cells from 
the freezer 



2 Grow cultures 
for desired 
number of 
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Option A: Grow two replicates of a 
senes of four 700 pi cultures from 
OD 0 0625 to OD 2 (mix replicate 
wells before harvesting) 
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Option B: Grow a series of two 50 ml 
cultures from OD 0 002 to OD 2 



3 Harvest 
cells 



4 Purify 
genomic DNA 



UPTAG DNTAG 



5 PCR amplify 
uptags and 
downtags 



6a. Hybridize to 
an array, wash, 
and scan 



3d 
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Figure 1. Workflow for pooled growth assay and barcode detection. Cultures are inoculated with thawed aliquots of pooled cells (step 1), and 
then grown for the desired number of generations (step 2) either robotically (Option A) or manually (Option B). Cells are harvested by 
centrifugation (step 3) and genomic DNA is then isolated from the harvested cells (step 4), uptags and downtags are independently amplified 
(step 5), and hybridized to an array (step 6a or sequenced directly step 6b). 



Copyright ©2011 Creative Commons Attribution License August 2011 | 54 | e2864 | Page 4 of 7 



IQve Journal of Visualized Experiments 



www.jove.com 



• NCP1 



ERG2 



* 



, ERG11-1 
ERG11-2 



° » 




7 



Genes (alphabetical) 



B 



in 



q . 




•BWP17 

• NCP1 

• ERG2 

• ERG11-2 

• ERG11-1 



0 5 10 15 20 
Time (hours) 

Figure 2. Sample data collected at certain points of the protocol. (A) Sample data from screening results (adapted from 13 ). The pool of tagged 
mutants was grown for 20 generations in the presence of clotrimazole and DMSO (control). Log2 ratio (control intensity/treatment intensity) was 
calculated and plotted as a function of gene. Highly sensitive strains (red) included the known target of clotrimazole, ERG11p. Note that this 
assay frequently uncovers other sensitive mutants in addition to the compound's actual target. Generally, these mutants are those that act 
synthetically with the target, those that are part of a general stress/treatment response, or are false positives that fail to confirm. (B) Example of 
confirmation data (adapted from 13 ). Results from the pooled growth assays can be validated by growing the strain in individual culture and 
compared against wild-type growth (black). 



common primers 




Illumina adapte 



DNTA6 region- ' • 

Illumina adapter 



Figure 3. Structure of the amplicon produced from pooled barcode assays for microarray hybridization or Barcode sequencing. The amplicon 
produced for each mutant in the collection contains homology to the genome for integration (blue regions labeled ATG and TAA), unique barcodes 
(labeled AG and indicated by a black dash). For microarray hybridization, the blue common primers are used to amplify a 60bp probe for 
microarray hybridization. For barcode sequencing, extended primers are used in the PCR reaction, comprised of sequences encoding the 
Illumina adapter (red bar), and index of 6bases cross-hatches) and the blue common primer for the upstream primer, and the same composite 
primer (minus the 6 base index) for the second primer. 
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Figure 4. Individual Growth assays for 1) Prescreening compounds against wild-type yeast to determine an appropriate dose for genome-wide 
screening and 2) Confirming results from genome wide screens. (A) A 96 well flat bottom plate is filled with 100 pi of cell suspension at an OD of 
0.062. Each well can contain the same strain (for dose-determination) or different strain and drug combinations (for confirmation assays). 2 pi of 
compound (typically dissolved in DMSO) is added and cells are grown with constant shaking for 16-20 h at 30°C. The final concentration of 
DMSO should not exceed 2%. In this example, in each well of the plate the growth curve is plotted in black against a plot of the control growth 
curve in red. (B) Higher resolution image of several prescreens obtained with an example drug overlaid on top of one another. In this titration 
series, an IC10-15 is obtained with the purple dose and would be appropriate for deletion profiling (HIP and HOP). Due to the non-linearity at higher 
optical densities, Tecan (or any similar plate reader) ODs must be calibrated using those obtained with a "traditional" 1mm path-length cuvette. 



Discussion 



Here, we outlined a protocol that, with modest modification, can easily be adapted to a wide range of existing collections of barcoded mutant 
collections of different microorganisms to create tagged mutant collections. We emphasize that while we have reported a protocol on the tagged 
transposon mutagenesis for the pathogenic yeast C. albicans, a very similar protocol could be adapted to a wide variety of unicellular fungi. 
Modified, this protocol works well in bacteria 13 , and currently collections for a number of additional fungal and bacterial genomes are under 
construction. At present, this assay provides the only comprehensive, genome-wide unbiased screen for gene-small molecule interactions. One 
particularly compelling feature of the assay is that no prior knowledge of the gene or small molecule is required. Despite the scope and power of 
these assays, their transferability to other labs has been hindered somewhat by the initial capital investment and informatics tools for the analysis 
of the results. With the adoption of a next generation sequence readout combined with robust tools for analysis, we expect their adoption to 
increase. 
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